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Abstract 

In this paper a unified approach for power allocation (PA) problem in multi-hop orthogonal frequency division 
multiplexing (OFDM) amplify-and-forward (AF) relaying systems has been developed. In the proposed approach, we 
consider short and long-term individual and total power constraints at the source and relays, and devise decentralized 
low complexity PA algorithms when wireless links are subject to channel path-loss and small-scale Rayleigh fading. 
In particular, aiming at improving the instantaneous rate of multi-hop transmission systems with AF relaying, we 
develop (i) a near-optimal iterative PA algorithm based on the exact analysis of the received SNR at the destination; 
(ii) a low complexity-suboptimal iterative PA algorithm based on an approximate expression at high-SNR regime; 
and (iii) a low complexity-non iterative PA scheme with limited performance loss. Since the PA problem in multi-hop 
systems is too complex to solve with known optimization solvers, in the proposed formulations, we adopted a two- 
stage approach, including a power distribution phase among distinct subcarriers, and a power allocation phase among 
different relays. The individual PA phases are then appropriately linked through an iterative method which tries to 
compensate the performance loss caused by the distinct two-stage approach. Simulation results show the superior 
performance of the proposed power allocation algorithms. 


I. Introduction 

Multi-hop relaying and orthogonal frequency division multiplexing (OFDM) are promising techniques for high¬ 
speed data communication among wireless devices that may not be within direct transmission range of each other. 
Relaying protocols are broadly categorized as amplify-and-forward (AF) relaying, in which each relay forwards a 
scaled version of the received noisy copy of the source signal, and decode-and-forward (DF) relaying, in which 
each relay forwards a regenerated version of the received noisy copy of the source signal. AF relays may also be 
categorized as blind/fixed gain, channel assisted, and channel noise assisted, based on how source-relay channel state 
information (CSI) and noise statistics affect the relay gains HI. The capacity analysis and transmission protocol 
design over relay channels have attracted lots of research activities in the past decade El- 0. As relays have 
shown their merit for data transfer purposes, multi-hop communications have been also included in advanced 
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wireless standards such as IEEE 802.1 In 13 , WiMAX, and LTE-Advanced 0 . To achieve power efficiency 

in multi-hop transmission systems, it is necessary to devise efficient power allocation (PA) strategies for the source 
and intermediate relays, when a multi-hop data transmission path is setup. For the simplest form of dual-hop 
relaying systems, the PA problem has been investigated in ifTOl - ifTSll . Specifically in IITSl . for a two-hop OFDM 
communication system with a given power budget, the authors presented the optimal power allocation at the relay 
(source) node for a given source (relay) power allocation scheme, that maximizes the instantaneous rate of the 
system. Then, using simulations they showed that by iterative power allocation between source and relay, a higher 
gain is achieved. A jointly optimal subcarrier pairing and power allocation scheme which maximizes the throughput 
of OFDM amplify-and-forward relaying systems subject to a statistical delay constraint is investigated in M- 
Power allocation for dual-hop OFDM relaying systems has been also considered in El- ED. The power allocation 
problem for relaying system models with more than two hops over narrowband fading channels has been considered 
in II 24 II . ||25]| . Especially in ESll . aiming at maximizing the instantaneous rate, the power allocation solution for AF 
relaying protocol over narrowband Rayleigh fading channels has been provided. In a path searching algorithm 
has been presented to find the best links among relays and then two subcarrier allocation algorithms are presented 
which aim at resource utilization improvement. Optimal and suboptimal power allocation schemes for multi-hop 
OFDM systems with DF relaying protocol has been developed in lIZTl . Adaptive power allocation algorithms for 
maximizing system capacity (when full CSI is available), and minimizing system outage probability (with limited 
CSI) have been proposed in ll28l for multi-hop DF transmission systems under a total power constraint. Aiming at 
maximizing the end-to-end average transmission rate under a long-term total power constraint, the authors in ll29l 
developed a resource allocation scheme for multi-hop OFDM DF relaying system in which the power allocated to 
each subcarrier and the transmission time per hop have been specified. In ||30l, the authors proposed the solution for 
power allocation problem in multi-hop OFDM relaying system (AF and DF) under total short-term power constraint, 
where the PA and capacity analysis is developed based on a high SNR approximation in the AF relaying protocol. 
The approximation used in BOl has been originally proposed in ll3Tl and performs well for small number of hops 
with high received SNR. In low-to-medium SNR regime or for multi-hop systems with more than three hops, this 
approximation loses its functionality in design of PA schemes. The multi-hop OFDM transmission system has been 
also considered in ||32l, where joint power allocation and subcarrier pairing solutions for AF and DF relaying 
protocols have been devised under short-term total power constraint. The analysis for AF relaying protocol in ll^ 
is also based on the high SNR approximation as discussed above. 

A. Paper’s contributions 

Although several research works have been reported on the power allocation problem in multi-hop OFDM systems, 
to the best of the author’s knowledge, no attempt has been made to develop a unifying approach addressing different 
aspect of these systems. In this paper, we focus on multi-hop OFDM systems with AF relaying and different system 
power constraints. We developed an unified framework for efficient power allocation which includes both iterative 
and non-iterative solutions. In particular, aiming at maximizing the instantaneous system rate under individual and 
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Source Destination 


Fig. 1. Multi-hop OFDM relaying system with A^-1 intermediate relays. 


total short and long-term power constraints, the following power allocation schemes have been devised: (i) A near- 
optimal iterative PA algorithm which is developed based on the analysis of an exact expression for the received 
SNR at the destination; (ii) A low complexity-suboptimal iterative PA algorithm in which we use an high-SNR 
approximation of the system rate for design purposes; and (iii) A low complexity-non iterative PA scheme based on 
a high-SNR rate analysis at the destination. The rest of this paper is organized as follows. In Section II we introduce 
the system model and the power allocation problem formulation. An iterative PA solution based on the analysis of 
the exact destination SNR is presented in Section III. In Section IV we provided sub-optimal PA schemes based on 
the analysis of an approximated high SNR expression at the destination. Simulation results are provided in Section 
IV, and concluding remarks are presented in Section V. 


II. System Model and Problem Formulation 

Fig. [T] shows the multi-hop transmission system model under consideration, where OFDM is utilized for broadband 
communication among the consecutive nodes. We assume that the system uses a routing algorithm, and therefore 
the path between the source and destination nodes is already established. Here the source node, Tq, sends data bits 
to the destination node T/v via A^-1 intermediate relay nodes, Ti, T 2 ,..., over orthogonal time slots and 

orthogonal subcarriers. The fading gain of the narrowband subchannel i (corresponding to the ith subcarrier) between 
nodes Tfc_i and T^, denoted by a^^i, is modeled as a zero mean circularly symmetric complex Gaussian random 
variable with variance In an AF multi-hop relaying system, each relay first amplifies the signal received from 
its immediate preceding node, and then forwards it to the next node in the subsequent time slot. The amplification 
gain in zth subcarrier of node is adapted based on the instantaneous fading amplitude of the channel between 
nodes Tk-i and Tk, i.e. \ak,i\- To ensure an output relay transmit power Pfc ^ on ith subcarrier, the amplification 
gain is adjusted as US: 


42 ^ Pk,i _ 


for fc = 1,..., A^ — 1; z = 1,..., Np 


where Pg.i ^nd Pj. ^ denote the transmission powers at source and fcth relay for zth subcarrier, respectively. In this 
model, the number of subcarriers, i.e., the number of points for fast Fourier transform (FFT), and the noise power 
at fcth relay within zth subcarrier are denoted by Np and Nq^. . respectively. 

For the considered multi-hop OFDM system model, the instantaneous received SNR over the zth subcarrier at 
the destination node is given by 

i't = (nL.d + *) -1) 


( 1 ) 
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where 7 ^ ^ denotes the instantaneous received SNR over the zth subcarrier of A:th hop with the 


average = 


I Using Q, the instantaneous rate of the end-to-end multi-hop system is given by 


t^=^Et"iiog(i + 77- 

> 

C may be rewritten as follows 


Considering ak,^ = and = 


^Ok 


c^ = FEfJiiog(i 


n£'=i(i+: 


)-l 


( 2 ) 


(3) 


PT<yk,-, 


In a Rayleigh fading environment, 7 ^ ^ follows an exponential distribution with the average ^ 

Given the power constraint, here the goal is to find the transmit powers of the subcarriers at the source and 
the relay nodes such that the instantaneous rate in (|^ is maximized. In this paper, we consider power allocation 
optimization problem under short-term power constraints (STPC), long-term individual power constraints (LTIPC), 
and long-term total power constraints (LTTPC). The general form of power allocation optimization (PA) problem 
in this work is as follows. 


max„ 


C 


s.t. 


E N {T N — 1 1 

i=l l^k=0 — f 


E(afc,i) = 


NxNf 


STPC 


LTIPC 


(4) 


= l LTTPC 

To explicitly identify the way power is distributed among different subcarriers and nodes, we denote by Pi and 
Pk,i, respectively, the allocated power to *th subcarrier (for all nodes) and the allocated power to Ah subcarrier in 
fcth node (hop). We also define two new nonnegative PA coefficients and f3k,i as follows 

Pi =fii X Pt (5) 

and Pk^i ^ ^ At ^ : ( 6 ) 

where ak,i = Pk,i x IJ-i- The optimization problem (3) may now be rewritten as follows with Hi and j3k,i as 
optimization variables, 

max/ 3 j^ C s.t. 

CL 

= l STPC 

E(At.) = 7, E(/3fc,,) = ^ LTIPC (7) 

ES E(aO = 1 Ef=”o' nPk,^) = 1 LTTPC 
C2. At — 0 z = 1,..., Np 
C3. /3fc,i>0 k = Q,...,N-l-i = l,...,NF. 

We note that the allocation of power to subcarriers (by finding /r^) and to subcarrier per node (by finding f3k,i) 
in problem ( 6 ) is equivalent to finding ^ in (3). Unfortunately, the power allocation problems (3) and ( 6 ) are 
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too complicated and cannot be solved with known optimization solvers. Hence, in the subsequent sections to hnd 
efficient solutions to 0 under different power constraints, we develop iterative and non-iterative algorithms based 
on the exact and approximate SNR expressions. 


III. Iterative Power Allocation 

In this Section, we consider the power allocation problem ( 6 ) as alternate maximization over two simpler power 
allocation optimization sub-problems. In the first sub-problem, the optimized {Pk,i} is obtained assuming that {/ii} 
is available. In the second sub-problem, {/i^} is determined for a given set of Next, we provide an iterative 

PA algorithm in which the two subproblems are alternately considered with the output of one as the input of the 
other. Numerical evaluation in Section V verihes the effectiveness of the proposed approach. 

A. Sub-problem 1: Power allocation among relays 

Given that the power allocated to each of the subcarriers {pi} is already known, the power allocation problem 
among relays in one subcarrier is equivalent to the power allocation problem in a multi-hop system with narrowband 
fading model. Here, we provide optimal solutions for the PA problems in multi-hop narrowband communication 
systems with individual long-term and total power constraints. We note that in an AF multi-hop transmission system, 
maximization of the instantaneous rate of the system is equivalent to maximizing the instantaneous received SNR 
1251 . The instantaneous received SNR at the destination of multi-hop system given in 0 may be expressed as 
follows: 


In (l + 7r') = + ^) = Eti ln(l + ^). 


Then, 7 ^ is found as 


7t = 


exp 


(EL.ln(l + i)) 


- 1 


1 -1 


( 8 ) 


Since, maximizing 7 ^ is equivalent to minimizing the argument of exponential function on the right hand side of 
we can conclude that 


argmax{^^ j 7 t = argmin{^^ j J2k=i ln(l + :^)- 
Using (|^, the optimization problem for the ith subcarrier under LTIPC is written as follows: 


(9) 


min{/3,„} Eto ln(l + 


( 10 ) 


s.t. E{l3k,i) = jj, k = 0,1,... ,N - 1, 

where {pi} is given. In appendix |0 it is shown that the objective function in ( [T0| ) is convex. Since the constraints 
are linear the main problem is convex, and therefore we use the Lagrange method to obtain the optimal solution. 
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The Lagrangian function is given by 

^N-l 


■C = ELo ln(l + + ELo (E(A) - *) 


( 11 ) 


where are Lagrange multipliers. By taking the derivative of C with respect to [3k,i, /c = 0,1,— 1, 

one obtains 


dC 

dh,i 


= 0 




l^k(k'ilk+i,iPk,i) 


= 0 


( 12 ) 


[3k,i — 


-l+,/l+4 




Substituting ( [T^ in the fcth individual power constraint in (10) leads to 

. -l+,/l+4^L2fc+E) 

E(/3fc,i) = ^ E 


r 

Jo 


1+4 


>^i1k + l,i 
^k,i 


2Aii7fc+i,i 




N 


(13) 


fk,rk+i,i ^7k + i,i - Al¬ 


to 

To find the Lagrange multiplier Ak,i, the integral in 0 can be numerically evaluated and a bisection root¬ 
finding method ll^ can be utilized. From ( fTS] ) and 0, it is evident that the power coefficient at the fcth node, 
[3k,i, k = 0,N — 1, is only dependent on the fading gain of its immediate forward channel. As a result, the 
proposed power allocation scheme can be implemented in a decentralized manner. Such a PA scheme is potentially 
applicable in ad-hoc wireless networks over narrowband channels. 

Using similar steps for the LTTPC case, the PA coefficient is obtained by 


-i+/+igE 

2/j.i'yk+i,i ’ 

where the constant A is calculated using the corresponding power constraint follows 

Af-l 

E = 1 - 

k=0 

Moreover, for the STPC case, the procedure is similar to the LTIPC scenario. 


(14) 


B. Sub-problem 2: Power allocation among subcarriers 

The instantaneous rate, C, can be written as a function of /iiS and f3k,iS, as follows 


nL7(i+7 


-)-i 


Given [3k,i^, C is to be maximized by finding the optimal piS. We start by expanding the expression for Ih 
follows 

n (1 + = 

k—Q 


(15) 


as 


N-l 


. , N-l N-l 

■ ) + if4(E 

N-l N-l N-l 

1/N\ / 1 _ 1 1 -I 

^ , n /3fc,i7)= + l,i /3j,i7j + l,i + 

k—0 3=0 1=0 

j:^k 


k-i'' Pk,ilk+1, 

k—Q 


( 16 ) 
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Then, one can rewrite the rate as 


C = 


^Et"iiog( 


1 + 




where 


N-l 

A _ J_ 1 

^l.i - 1! 2^ I3k,.lk+i,: 

k—0 

N-l N-l 

^ 2 ,i = ^ S S 


2! ^ ^ 8k i7fcj.i - 

fc =0 i -0 + 


1 

ft-.i7j + l,i 


(17) 


( 18 ) 


^ N-l N-l N-l ^ ^ ^ 

^N,i — jv! S S 0k,ilk+i,i /3i,i77 + i.i /3(,i7i + i,i ■ 

Under LTIPC, we construct the following Lagrangian function 

- A*(E(^j) - ]^)) - V^fl, 

where v’^ = [ui,..., and /x = [^i, ■ ■ ■ and Ai, and Vi are the Lagrange multipliers corresponding 

to the constraints Cl and C2 in 0, respectively. Furthermore, by taking the derivative of C with respect to 
fii, i = 1 ,..., Np, one obtains 

(A, + u,)D^ + ((A, + d' = 0, (20) 

where 

D = ^ ^ + ^Ar,i 

and D = dD/d^i. 

To find the PA coefficient the polynomial equation in ( |20| ) is to be solved numerically, by considering the power 
constraint. The difficulty for solving ( (20l i increases when the number of hops increases. For STPC and LTTPC, the 
procedure is the same as LTIPC, with their corresponding power constraints. 


C. The iterative scheme 

In this subsection, we present an algorithm which iterates between power allocation among relays and sub¬ 
carriers to maximize the instantaneous transmission rate of the network. As verified in Section VI, such an 
iterative solution can be used as an upper bound for the performance evaluation of the proposed suboptimal 
solutions in Section IV. In this algorithm, optimizations in sub-problem 1 and sub-problem 2 are repeated al¬ 
ternately, such that the PA parameters obtained from the previous optimization are the input for the next one. 
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Iterative Algorithm 

1. Initialize the subcarrier PA coefficients to Hi = 
l/NFyi = l,...,NF. 

2. Given ^^s, hnd (3k ,from the sub-problem 1. 

3. Find the instantaneous rate C by substituting the PA 
coefficients (3k, i and in ( [TS] ). 

4. Given (3k,i&, find /i^s from the sub-problem 2. 

5. Find the instantaneous rate C in ( [T5] l using (3k,i and 

obtained in steps 2 and 4, respectively. 

6. If the difference between rates found in steps 3 and 
6 is above a given small value, repeat steps 2 to 6, 
if not (or after a predefined number of iterations) go 
to step 7. 

7. Report and (3k,i, Vi = 1,, Nf ,fc = 1,... , N — 

1 . 


IV. Power Allocation Schemes in High-SNR Regime 
I n this Section, we focus on a high SNR regime in the multi-hop system and present efficient solutions for 
the power allocation problem in (6). The motivation to study such power allocation algorithms is their simplicity 
with respect to the iterative solution presented in the previous section. We start by rewriting the instantaneous rate 
expression in GD as 

c^ = ^EfJiiog 

In high SNR regime, the parameters Ak,i in (17), for fc = 2, • • • ,N, are negligible in comparison with Ai,i, hence 
we neglect higher order terms and rewrite the expression in ( |2T| l as follows 

Cap, = jrJ:l\log[l + ^). ( 22 ) 

Substituting Ai,i from ( |T8] l, the instantaneous rate is expressed as 

Capp = Efci log(l + ^ -)■ 

Using (|2^, iterative and non-iterative power allocation schemes are developed in the following subsections. 



( 23 ) 
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A. Iterative PA in high-SNR regime 


In high-SNR regime, the steps 2 and 4 in the iterative algorithm can be implemented in a simpler way, as 
described below. 

A.l Sub-problem 1: power allocation among relays 

According to (|2^, the instantaneous rate of the ith subcarrier is given by 


Ciapp — log(l -f ^ ^ ). (24) 

(So 

By formulating the power allocation problem for multi-hop narrowband system, the general optimization problem 
is written as follows: 

app 

S.t. 

ESTo' = 1 STPC (25) 

E(/3fe,,) = S( fc = 0,l,...,A^-l LTIPC 

ES:o'E(/3m) = 1 LTTPC 


which is a convex power allocation problem and using the Lagrange method, the PA coefficients under long-term 
individual power constraint are derived as 


l^k,i — 


N 


1 _ k(=+i.d 


And, the power allocation coefficients under LTTPC are calculated as 

/3fc,i = ^ ^ 


Af E 

j=i 


^ |“(=+i,i| 

,S^7^ 


(26) 


(27) 


Moreover, the power allocation coefficients under STPC are derived as follows 

(28) 

Note that the PA solutions in (|26]l-(|28]l were first derived in ll25l . where the authors investigated high-SNR power 
allocation problems for a multi-hop narrowband system model. However, here we use such solutions to provide 
power allocation scheme for a wideband multi-hop system with OFDM modulation. 

A.l Sub-problem 2: power allocation among subcarriers 
In high-SNR regime, the optimization problem in Q is rewritten as 



Capp S.t. 

Cl. 

E.Ei = 1 STPC 

E(/i,) = ^ LTIPC 

EtlE(Ai^) = l LTTPC 

C2. Pi ^ 0 i = 1,..., Np. 


( 29 ) 
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where Capp is given in ( |2^ . Since the objective function in ( [29l l is concave (see Appendix [B]l and the constraints 
are linear, the optimization problem (23) has a unique optimal solution. 

Under LTIPC, one can construct the Lagrangian function as follows 


^ = EfJiiog(i 


E N — 1 
k~' 


Mi 




-) 




(30) 




where A^, i = 1,2,Np, and the vector = [rii,..., are Lagrange multipliers. The optimized is 


calculated by setting = 0, as follows 


Mi = 


N-1 


1 _“ 


Ai /l^fc=0 7fc+i,i/3fc,, 

In (25), the constant Ai is found to satisfy the constraint Cl in ([29ll with equality, i.e.. 


E 


N-l 


1 _ 

Ai 2—ik=' 


or 


where 


fe =0 7fe+i,i/3fc,i 


fY{y)dy = 


N-l 


1 

N fr 




(31) 


(32) 


(33) 


(34) 


■'k=0 'ik+i,iPk,i ■ 

In ( |33| ), fviy) is the probability density function (PDF) of Y. The PDF of Y is calculated by the convolution of 
probability density functions for , k G {!,■ ■ ■ ,A^—1 }, where 7 fe+i,i follows an exponential distribution, 

and /3fc i is a given constant. 


Using the same procedure, the solution for this sub-problem 2, under the LTTPC is derived as 

n -I- 


Mi = 


tv-i 


A 2^k— 


A ^k=0 'ik + i,iPk,i 

where the constant A is found to satisfy the constraint Cl in (|29ll with equality, i.e.. 


EEf 


1 + 


7)! + l,i/3fc,i 


= 1 . 


(35) 


(36) 


The procedure for STPC is the same as LTTPC. 


B. Non-iterative PA scheme in high-SNR regime using channel statistics 

However in previous sections we have presented two PA schemes which utilize CSI for power allocation among 
subcarriers and relays, low complexity schemes which can work with channel statistics instead of the CSI are always 
appreciated. As we saw in Section |IV-A[ specially in (|26ll-(|28ll, power allocation among relays is independent from 
power allocation among subcarriers. This fact motivates us to present a non-iterative power allocation algorithm in 
high-SNR regime. To this end, we insert the derived in ([26]l-(|28]l into (|3T]i and ([35]l. Then, the PA coefficient 
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C(k,i = IJ‘iPk,i is derived by substituting the instantaneous values with their means. In this case, the power allocation 
coefficient under long-term total power constraint will be as: 

N-l r- N 

A “ 


^k,i 


N-l 


n + r 


TT 1 

^ . . Y j- 


N 

E 


n -1 




r • • 


/c = 0 V 

For the case of long-term individual power constraint, the solution is given by 

\/r k+l,i 




N-l 




Ny/TT 


— - Yl 

V7fc+i.irfc+i,i 


Ny/TTJk+l,t ' 


Moreover, by considering a short-term power constraint, we obtain 


^kA — 


■ A-l N 

V V 

1 

+ 

N 

/ t'o VT'fc+i.* ^ 

/7i,i 


1_ 


1 -1 


where A and A; are found using STPC, LTTPC, and LTIPC power constraints in 
section we evaluate the performance of proposed power allocation algorithms. 


respectively. In the next 


V. Performance Evaluation 

This Section presents simulation results for performance evaluation of the proposed power allocation schemes. In 
simulations, we assume a multi-hop AF relaying system over Rayleigh fading channels under short and long-term 
individual and total power constraints. In all following hgures, EPA, ASY, IT-EXA, and IT-ASY refer, respectively, 
to the equal PA method, the PA scheme using high-SNR analysis (Section IV), the iterative power allocation 
algorithm using exact SNR analysis (Section III), and the iterative power allocation algorithm using high-SNR 
analysis (Section IV). 

Fig. 1^ shows the average rate of the 2-hop OFDM relaying system with Np =64 subcarriers versus average SNR 
of the direct link Fq, under long-term total power constraints. Balanced and unbalanced links are considered here. 
For modeling unbalanced links in multi-hop system, we adopt the setup of ll25l . in which it is assumed that the 
/cth terminal, k = 1, ...,N + 1, is located at the distance dk = from its previous terminal, where d 

is the distance between source and destination. Hence, using the Friis propagation formula lf34l . the average SNR 
of /cth hop is given by F^ = ( (zv-fF^a''+ 2 ) where 5 is the path loss exponent and Fq is the average SNR of 
direct link. We consider 5 = 4 in this work. For balanced links, the inter-distance among nodes is considered the 
same, then the average SNR of fcth hop is given by F^ = (TV + l)^Fo. 

In Fig. 1^ one can see that the iterative scheme (IT-EXA) acts as an upper-bound for other power allocation 
schemes. The iterative scheme using high-SNR analysis (IT-ASY) has acceptable performance with a lower com¬ 
putational complexity because we have closed-form expressions for fips and iipS in IT-ASY scheme. Also the 
non-iterative scheme using high-SNR analysis (ASY) shows a superior performance with respect to equal power 
allocation, however, the difference between ASY and IT-EXA scheme is quite large because we use channel statistics 
instead of the CSI in ASY scheme. 

Fig. 0 shows the average rate of the 3-hop OFDM relaying system with 64 subcarriers versus average SNR of 
the direct link, Fq, under short-term power constraints. One can see that the performance of iterative algorithm 
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based on exact SNR analysis is superior than other schemes. Moreover, the performance of the non-iterative power 
allocation scheme is close to that of the low complexity iterative algorithm, and both of these schemes provide 
significant gain over the equal power allocation. The same results are depicted in Fig. for 2-hop system under 
long-term total power constraint. 

In Fig. 1^ the outage probability of 2-hop OFDM system (Np = 64) is depicted versus the average SNR of 
direct link, Fq, for several power allocation schemes under long-term individual power constraints. The outage 
probability is defined as the probability that the instantaneous rate falls below 1 bit/sec/Hz. In Fig. the average 
rate of 2-hop OFDM system when iterative PA algorithms in both exact and asymptotic forms are applied (under 
LTTPC) is depicted. In this figure, the index i stands for the number of iterations considered, while i = 0 refers 
to the uniform power allocation scenario. As an interesting observation, from this figure we can see that after a 
few iterations the average rate converges to its maximum value. We also note that from a complexity perspective, 
the computational complexity of proposed iterative algorithms is directly related to the complexity of sub-problems 
in each iteration. As an example, for the 2-hop OFDM system considered in Fig. the complexity of IT-ASY 
and IT-EXA algorithms may be easily related to the complexity of solving PA sub-problems among subchannels 
and relays (the corresponding waterfilling solutions for these sub-problems are presented in Section III and IV). In 
particular, the complexity of waterfilling solutions has been already investigated in 051 . 

From Fig. |2]|^ we make the following observations: (i) The IT-EXA scheme provides the best performance 
among other methods at the cost of a higher computational complexity; (ii) The IT-ASY scheme provides a data 
rate performance which is very close to that of the IT-EXA scheme, while it enjoys a considerable lower complexity; 
(iii) The ASY scheme provides an acceptable level of data rate performance due to its considerable lower complexity 
and easier implementation as it needs channel statistics instead of CSI; (iv) The performance of the IT-ASY scheme 
in high SNR regime converges to that of the IT-EXA scheme, as it is verified in Fig. and Fig|^ (iv) For a 
multi-hop OFDM scenario, the proposed power allocation schemes greatly outperform the scheme with uniform 
power allocation. 


VI. Conclusion 

In this paper we considered the problem of power allocation in narrowband and broadband (OFDM) multi¬ 
hop relaying systems employing non-regenerative relays with different power constraints. We proposed exact and 
approximate design approaches depend on the wireless application demand and network structure. In particular, 
aiming at maximizing the instantaneous multi-hop transmission rate, several power allocation algorithms have been 
developed in a unified framework including: (i) an iterative power allocation method which provides an upper- 
bound performance; (ii) a relatively low-complexity iterative power allocation method; and (iii) a non-iterative 
power allocation scheme with acceptable performance at high SNR regime. Moreover, we provided performance 
comparison with respect to an equal power PA solution and quantify the rate performance loss incurred at the price 
of low complexity and low feedback overhead. 
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Fig. 2. Average rate of 2-hop OFDM {Np = 64) relaying system under LTIPC. 


Appendix A 


Here, we prove that the objective function in (fTO)! is convex. Let denote this function by f(/3fe) as follows 


TO = Er=YMi + 


/3fc7fc4 


We can easily obtain 


d’^fjPk) _ l-F2/3fc7fc + i 

(/3fc(/3fc7fc+i + l))^ ■ 


(37) 


(38) 


Since the coefficient jSk is between 0 and 1 (see Section II), and 7^+1 takes positive values, the second derivative 
in is positive for any channel realization, and fiPk) is convex. 


Appendix B 

Here, we show that the objective function in ( |29| ) is concave. Let f(pi) denote the objective function, that is 

I 1 I ^jPT 


f(A*) = Efc=i log 1 1 + 777777;)- 




(39) 


We rewrite this function as follows 

f(Ai) = Er=i log (1 + AiffO 


( 40 ) 
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Fig. 3. Average rate of 3-hop OFDM {Np = 64) relaying system under STPC. 


where 


Si = 


Pt 




'‘Ok.i 


Taking the second derivative with respect to I3k, one obtains 


(41) 


log(e)gi^ 

As stated in Section II, ni is between 0 and 1 for any channel realization. Hence, the second derivative in ( |42] l is 
negative for any channel realization, and the function f(/ii) is concave. 
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